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Wave Solutions 



Overview 

In classical continuum physics, a wave is a mechanical disturbance. Whether the dis- 
turbance is stationary or traveling and whether it is caused by the motion of atoms 
and molecules or the vibration of a lattice structure, a wave can be understood as a 
specific type of solution of an appropriate mathematical equation modeling the un- 
derlying physics. Typical models consist of partial differential equations that exhibit 
certain general properties, e.g., hyperbolicity. This, in turn, leads to the possibility of 
wave solutions. Various analytical techniques (integral transforms, complex variables, 
reduction to ordinary differential equations, etc.) are available to find wave solutions 
of linear partial differential equations. Furthermore, linear hyperbolic equations with 
higher-order derivatives provide the mathematical underpinning of the phenomenon of 
dispersion, i.e., the dependence of a wave's phase speed on its wavenumber. For sys- 
tems of nonlinear first-order hyperbolic equations, there also exists a general theory for 
finding wave solutions. In addition, nonlinear parabolic partial differential equations 
are sometimes said to posses wave solutions, though they lack hyperbolicity, because 
it may be possible to find solutions that translate in space with time. Unfortunately, 
an all-encompassing methodology for solution of partial differential equations with any 
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possible combination of nonlinearities does not exist. Thus, nonlinear wave solutions 
must be sought on a case-by-case basis depending on the governing equation. 



Hyperbolic Equations 



l|. It is undeniable, however, that we 



A wave solution is difficult to define precisely 
know a wave when we see it. Mathematically, partial differential equations (PDEs) of 
hyperbolic type have certain properties that allow them to possess wave solutions in 
a rigorous sense. For convenience, henceforth, we assume all equations are properly 
nondimensionalized so that we need not speak of units and dimensions. We denote by 
x = (xi, . . . ,XN d ) the spatial coordinates in ^-dimensional Euclidean space and by t 
the temporal variable. 

Balance Laws 

Consider the change in some field quantity u = u(x, t) (e.g., displacement, velocity, 
electric or magnetic field, etc.) over a domain Q. When u is conserved, the rate of 
change of its volumetric average must balance the flux through (in and out of) the 
boundary dfl of the domain and also the production/loss of u within the domain: 



M V L Ud " = " V L g A ' F ' <U) ^ + V L R(U) * 



(1) 



where dv is an infinitesimal volume element, V is the volume of i?, n = (ni, . . . , fiN d ) is 
the surface normal to df2 and da its infinitesimal surface element. Equation ([[]) is called 
the integral form of the balance law for u with {Fi, . . . , F^} and R, respectively, being 
the fluxes in the coordinate directions and the source term. 

Assuming, for the sake of argument, that all mathematical prerequisites are 
satisfied, we can make use of the divergence theorem and require that Eq. ([T|) holds for 
any fixed (stationary) domain Q, to obtain the differential form of the balance law: 
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When R = 0, we say that Eq. (j2J) is a conservation law. Furthermore, this system is 
hyperbolic if the composite Jacobian J(u; n) = Yli^j %"^T"( U ) nas om y rea ^ eigenvalues 
and is diagonaliazable for all n and u 2|, Chap. 18]. If the eigenvalues are also distinct, 
then the system is called strictly hyperbolic. The eigenvalues A; and right-eigenvector 
may both depend on u as well; thus, if • ^ for all u, the Zth eigenpair is 
said to be genuinely nonlinear, while it is linearly degenerate when the latter condition 
fails to hold. It is important to note that the definition of hyperbolicity rests upon the 
assumption that each Fj does not depend on Vu (or higher-order gradients) explicitly 
because such a dependence would not allow us to define a Jacobian as done above. In 
other words, this notion of hyperbolicity only applies PDEs that can be expressed as 
first-order systems. 

Linearization 

By introducing dimensionless groups of parameters, some of which are assumed to be 
small, one can formally develop a solution of a nonlinear equation into an asymptotic 
expansion. At the leading order, the equation to be solved is linear and independent of 
said small parameter (s). Alternatively, the problem can be "linearized" by evaluating 
all nonlinear terms at some equilibrium state, which is usually equivalent to the formal 
asymptotic expansion. Linearizing Eq. (j2J) leads to 

j=l 3 

where Aj := dFj/du and B := dH/du represent Jacobian matrices evaluated at some 
(constant) equilibrium state u . For convenience, we have assumed that R depends 
only on u. Moreover, u has been redefined to eliminate all constant terms on both 
sides of Eq. ©. 
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One-way Wave Equation 

In the one-dimensional case (N d = 1), leaving the j = 1 subscript understood, the sys- 
tem in Eq. fl3]) can be decoupled whenever A and B are simultaneously diagonalizable, 
i.e., A = S _1 AS and B = S~ 1 TS with A and r being diagonal matrices with the 
eigenvalues of A and B along their diagonals. Then, Eq. ([3]) is equivalent to 



where u = Su, and / ranges over the number of components of u. Because Eq. (j3J) has 
the same form for any I, we drop the I subscript; tildes are also left understood. 

It is now a simple exercise to show that Eq. (j2J), subject to u(x,t = 0) = u(x), 
has the following progressive wave solution: 



The phase velocity c of the wave is determined by A, while its attenuation is determined 
by 7. For A > (resp. A < 0), u(x,t) is a right (resp. left) traveling wave. As a result, 
Eq. (j3J) is often referred to as the one-way wave equation with a source term. When 
7 = (i.e., there is no source in the original balance law), the initial wave form 
propagates unchanged. Notice that when 7 > there is unbalanced "production" of u 
in the domain, leading to the unbound growth of the solution in time. When 7 < 0, on 
the other hand, there is unbalanced "loss" of u within the domain, leading to u(x, t) — > 
for all x as t — > 00. The solution given by Eq. (jSj) is illustrated in Fig. QJa). 

An important feature of the solution given by Eq. (j^J) is its finite domains of 
influence and dependence. For any given initial condition, the value u(xq) for all Xq 
in the domain can travel only as far as x = Xq + ct during the time interval [0,t], 
so its domain of influence is the spatial interval [xq,x + ct], where we have assumed 
c > without loss of generality. Similarly, the value of any solution u(x,t) for all x in 



at ox 



dui dui 



(4) 



u(x, t) = e 7 *-u(x — ct), c = A 



(5) 
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the domain can only be influenced by those in the spatial interval [x — ct,x], which is 
called its domain of dependence. Both of these domains are intervals of finite length ct. 
Thus, Eq. (j3J) exhibits finite speed of propagation of signals at a phase velocity given 
by c. Since all hyperbolic systems of first-order hyperbolic systems can be linearized 
and written in the form of Eq. Q, which can be decoupled as in Eq. we say that 
hyperbolic PDEs posses wave solutions in a strict sense. 

Two-way Wave Equation 

Again in the one-dimensional case, suppose u has two components, while A has only 
off-diagonal entries. Leaving the tildes understood, Eq. (j3J) is then equivalent to the 
following coupled system: 

+ A u ^r = Bxxu x + B 12 u 2 (6a) 
at ox 

^jr + ^ 21 ~^T = B 2i«i + B 22 u 2 (6b) 

Furthermore, suppose that B = 0, then u 2 can be eliminated from the above system 
(by taking the t derivative of the first equation and the x derivative of the second 
equation) to obtain a scalar equation for u = U\. 

cPu = 2 cPu 

dt 2 dx 2 y } 



where c := \J A\ 2 A 2 \. Note that A\ 2 A 2 \ > by the assumption that the system in 
Eq. © is hyperbolic. It can be shown that Eq. (jJJ), subject to u(x,t = 0) = u(x) and 
^(x,t = 0) = v(x), has the following progressive wave solution: 

1 1 r x + ct 

u(x, t) = - [u(x — ct) + u(x + ct)} H / v(s) ds (8) 

2 2c J x _ ct 

This solution of the initial-value problem on — oo < x < +oo, t > is the celebrated 



d'Alembert solution of the wave equation 



^4-1]. Note that, for any given set of initial 



conditions {u, v}, the domain of influence at any xq is [xq — ct, Xq + ct]. Meanwhile the 
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domain of dependence of a solution at any x is [x — ct, x + ct]. Both of these are finite 
intervals of length let. Since the two-way wave equation allows propagation of signals 
in two directions, the domains of influence and dependence are twice the size of those 
of the one-way wave equation. The solution given by Eq. (jHJ) is illustrated in Fig. [Hh>). 

Damped Wave Equation 

Returning to Eq. fl6]), let us now take B\\ ^ while Byi = -B21 = -B22 = 0. Eliminating 
u 2 between the two equations now leads to 

d 2 u Idu 2 d 2u 

'W + r'dt =C dx^ O 



where c := y/A^A^i and r := —1/Bu. The term proportional to r, which is usually 
called the relaxation time, corresponds to damping (resp. growth) when r > (resp. 
r < 0). Equation is also known as the telegrapher's equation. It arises in heat 
conduction when the heat flux has exponentially-decaying "memory" of the history of 
the temperature gradient |4|, [5|. Equation (jHJ) possesses a d'Alembert-type progressive 
wave solution [3|, §5-5]: 

u(x, t) = - [u(x - ct) + u(x + ct)} e ~ t/(2r) 



where p(s) = p(s; x, t) :— y 1 — (x — s) 2 / (ct) 2 and /„(•) is the modified Bessel function 
of the first kind of order v. Clearly, the damped wave equation has the same domains of 
influence and dependence as the two-way wave equation. The solution given by Eq. (TTOT) 
is illustrated in Fig. [TJc). 
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Fig. 1. Progressive wave solutions u(x,t) (superimposed at four consecutive times), starting with a 
Gaussian initial condition u (dashed curves) and v = 0, of the (a) one-way wave equation (7 = 0), (b) 
two-way wave equation and (c) damped wave equation (t = 5c, dotted curves show the exponentially- 
decaying amplitude of the peaks). 

Standing Wave Solutions 

The two-way wave equation on a finite domain, unlike the one-way wave equation, 
allows for the possibility of waves to superimpose in a way such that a standing wave 
pattern forms. To illustrate this phenomenon, consider the two-dimensional two-way 
wave equation on a domain Q: 

d 2 u 

— = c 2 Au, x e n, t>0 (11) 
dt 2 v ' 

where A := V • V is the Laplacian operator. Furthermore, due to the presence of a 
second-order spatial derivative, we must specify the value of u(pc,t) everywhere on the 
boundary df2. 

For such initial-boundary-value problems, the technique of separation of vari- 
ables |3|, Chap. 7] can be used to seek a solution of the form w(x, t) = X(x.)T(t), 
where X and T are to be determined. This assumption reduces the problem of stand- 
ing waves to one of finding the appropriate set of eigenf unctions {A' m (x)}, where 
m = (m 1; . . . , mjy d ) is an integer multi-index, that satisfy — A^ m (x) = yU m A' m (x) 
in Q and the appropriate boundary conditions on df2. Using linearity of the PDE 
and the definition of eigenfunctions, the task of finding T(t) is reduced to solving a 
second-order ordinary differential equation (ODE) with constant coefficients. Then, the 
general solution to Eq. f TTTj) is a superposition of all eigenfunctions: 
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x, t) = ^ [a m sin(yTi^ct) + /3 m cos(yT^ct)] # m (x) (12) 



where the coefficients {a m , /3 m } must be determined from the initial conditions w(x, 0) 
and ^r(x, 0) using orthogonality relations for the eigenf unctions. Unlike the progressive 
wave solutions above, the solution in Eq. ([12]) does not depend on the wave coordinate 
x — ct, rather it is standing wave — a superposition of purely spatial patterns, given by 
the eigenf unctions {A'm( x )}) whose coefficients are evolving in time. 

A thorough treatise on eigenfunctions, orthogonal systems and series expansions 



is given by Courant and Hilbert 
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Chap. II]. See also the discussion of Eq. (lTTj) in 
Chap. V §5], where it is shown that all /i m > for the Laplacian operator, justifying 
the form of the solution for T m {t) used in Eq. f|T2|) . 

Dispersion 

For linear (not necessarily hyperbolic) PDEs with constant coefficients posed on an in- 
finite domain, one can seek Fourier- mode wave solutions of the form w(x, t) = e * k x - tw *. 
Here, we only consider the scalar case, however, the approach is readily extended to 
systems of linear equations [7|. For k fixed, e lk ' x_;iaj ' is called a monochromatic wave 
solution. Upon substitution into the governing PDE, the problem is reduced to solving 
an algebraic equation of the form S(oj, k) = 0, which is called the dispersion relation. 
The solutions for u as a function of k, one or more of which may be complex-valued, 
determine the allowed wave modes of the PDE. 

As an example, consider the following one-dimensional linear PDE with constant 
real coefficients a, b and c: 

du du ^d 2 u d 3 u 
dt dx dx 2 dx 3 

The resulting dispersion relation is 

E(w, k) = -iu + aik + b(ik) 2 + c{ikf = (14) 
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Solving for u>, we immediately obtain u(k) = ak + ibk 2 — ck 3 , which can be interpreted 
as follows. For a solution of the form e * k ' x_;iaJ ', the phase velocity (sometimes denoted 
as c p ) has the natural definition c := 3f?[w]k/|k| 2 = ^R.[u]/k for one- dimensional propa- 
gation. Clearly, the third-order spatial derivative in Eq. ( TT3|) leads to a c p that is not 
independent of k, i.e., modes with different wave numbers have different frequencies 
(and, as a result, different phase velocities). This is the phenomenon of dispersion. 
Meanwhile, the second-order spatial derivative in Eq. (1131) leads to u having an imagi- 
nary part, which results in either growth (b > 0) or decay (b < 0) of the mode because 
e ~iut = e 9Mt e -tRMt_ If gq^j = o (e.g., a = c = in Eq. ([13])), then the PDE does not 
possess wave solutions, strictly speaking. 

To fully describe dispersive phenomena, the concept of group velocity c g : = 
Vk9?[w(k)] must also be introduced. For one-dimensional propagation, we simply have 
c g = d$t[(J\/dk. To understand the difference between phase and group velocities, con- 
sider a superposition of two one- dimensional monochromatic waves with nearby wave 
numbers given by k and k + Ak, where Ak <C 1. Then, in a dispersive medium, we 
have 



( x ^ _ e i[kx-w(k)t] _j_ e i[{k+Ak)x-ui{k+Ak)t] ^ e %k{x-[ui{k)/k]t} ^ _|_ e iAk{x-[duj/dk]t}^ 



where we have used the truncated Taylor expansion co(k + Ak) w w(k) + ^Ak and 
taken u(k) to be real without loss of generality. The first term on the right-hand side of 
Eq. (TT5l) represents a wave traveling with phase velocity c p = u(k)/k. However, due to 
dispersion (i.e., the fact that u(k)/k is not constant or, equivalently, du/dk ^ 0), the 
wave now has a modulated envelope that travels with its own phase velocity, namely, 
the group velocity c g = %. Thus, we arrive at the following intuitive definitions: the 
phase velocity describes the coherent propagation of a single mode with wavenumber 
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k, while the group velocity describes the coherent propagation of a group of modes 
with nearby wavenumbers in the range [k — Ak, k + Ak]. 

In general, if we were to project the initial condition onto the Fourier modes, 
use the dispersion relation to eliminate u in favor of k and "sum" over the continuum 
of wave numbers (i.e., integrate over k), then we would obtain the general solution to 
the linear PDE. What we have just described is the Fourier transform as a solution 
technique on an unbounded domain. A thorough discussion of linear dispersive wave 
phenomena is given by Whitham 8|, Chap. 11]. 

An important application of dispersion relations is in the study of instability. 
When a nonlinear PDE is linearized the coefficients carry some meaning about an 
equilibrium state (as discussed above). If the dispersion relation allows for solutions for 
uj with $s[ou] > for certain combinations of parameters (e.g., b > for Eq. f JT3|) ). then 
exponential growth of a Fourier mode can occur. Seeking to determine whether such 
modes exist is the subject of linear stability analysis, which has important implications 



in the theory of thermal stresses 



9j, Chap. 9]. In practice, such an analysis must also 



take into account specific boundary conditions on finite domains. 
Nonlinear Solution Types 

For simplicity of presentation, let us restrict to one-dimensional scalar conservation 
laws: 

| + s F <"> = ° (16) 



All equations of this form, subject to some 



formally) by the method of characteristics 



given 
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initial data, can be solved (at least 



Method of Characteristics 

Let us first rewrite Eq. ffT6]) in quasilinear form: 
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i + F'(u)4-}u- (17) 



dt dx 

The physical interpretation of Eq. ( ITT)) is that w is advected by the "velocity field" F'{u) 
and, since this "material derivative" vanishes, u does not change with time in the co- 
moving frame. Therefore, for a given initial condition u(xq,0) = u(xq), u(x,t) = u(xo) 
for all x and t such that x = X(t), where X(t) is the solution of the initial- value 
problem 

— = F'(u(X(t),t)), X(0)=x (18) 
The curves X(t) are a called the characteristics of Eq. (JT6J) - The analysis is easily 
extended to systems of hyperbolic conservation laws by using the fact that the Jacobian 
dF/du is assumed diagonalizable, thus the system of equations can be decoupled locally 
by a change of dependent variables, as we showed earlier for linear PDEs. 

Since u is conserved along the characteristic curves, we easily see that the general 
solution of Eq. (fTST) is 

X(t) = x + F'(u(x ))t (19) 
To evaluate u(x, t) we have to solve the (possibly) nonlinear equation x = X(t) for x ; 
i.e., we have to trace back along the characteristic. Notice that the characteristics given 
by Eq. (TT9|) are, in fact, lines because F' does not explicitly depend on x or t. 

Depending on the functional form of F', different characteristic curves can have 
different slopes. For simplicity, consider the case of u being a step function. This initial- 
value problem for Eq. (TT6l) is termed the Riemann problem. Then, depending on the 
left and right values u T := lim^o? u(x) of this step-function initial condition, three 
distinct types of characteristics exist as shown in Fig. [2] [see also 8, Chap. 2]. 

First, in Fig. |2(a), if F'(u + ) > F'(u~), characteristics with slopes that continu- 
ously vary between F'(u + ) and F'(u~) must emerge, which gives rise to a rarefaction 
fan. It can be shown that u(x,t) must vary continuously between u + and u~ in the 
fan. Second, in Fig. [2](b), if F'(u + ) < F'(u~), then the characteristics carrying these 
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two values cross. Beyond such crossings, the solution u(x,t) is multivalued. To resolve 
this difficulty, a shock wave is introduced, i.e., a solution that transitions discontin- 
uously between the values u T where the characteristics cross. Third, in Fig. [2Jc), if 
F'(u + ) = F'(u~) but u + 7^ u~ , then the characteristics remain parallel but the one 
emanating from xq = separates the two (initial) values of u for all time. Such a wave is 
called a contact discontinuity and is encountered in hyperbolic systems of conservation 
laws for which an eigenpair is linearly degenerate. 




t 



Fig. 2. Characteristics in the (x,t) plane corresponding to a Riemann problem initial condition for 
(a) a rarefaction wave solution (leading and trailing edges of the rarefaction fan are bold), (b) a 
shock wave solution (crossing characteristics are replaced by bold shock position) and (c) a contact 
discontinuity wave solution (separating characteristic is bold). 

Rarefaction Waves 

Since all initial condition emerging from a step-function have Xq = 0, then we deduce 
from Eq. (fT9j) that the characteristics in a rarefaction fan satisfy x = F'(u(x,t))t, or 
F'(u(x, t)) = x/t. Since F' is a function of u alone, it follows that u must be a function 
of the similarity variable x/t. In other words, rarefaction waves are self-similar solutions 
of Eq. f|T6|) in terms of x/t. 

Shock Waves 

When the characteristics cross each other, the solution jumps discontinuously between 
two values. Hence, it does not posses finite derivatives across the shock, and cannot 
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satisfy Eq. (fl6l) as written. Thus, more generally, the concept of a weak solution u(x, t) 
of Eq. (j2J) can be introduced. Weak solutions satisfy 



I J n [i + £^F j( u) -R(u) j^i)^* = 



(20) 



for all "smooth enough" (p(pc,t) that vanish identically outside the finite space-time 
domain Q x (0,T). All strong solutions, i.e., those that satisfy Eq. 02]) directly, are 
also weak solutions. However, weak solutions are not unique as there can be many that 
satisfy Eq. (120]) and the given initial/boundary conditions. 



To determine the unique weak solution, one appeals to the Lax condition 



For a system such as Eq. (Tl6|) . this condition states 
physical if there exists an I such that 
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12|. 



Chap. 10 §3] that a shock is 



A,(u-) > V > A,(u+) (21) 

while 

Vi(u-) < V < A m (u+) (22) 
In other words, the Lax condition states that a shock is physical only if "information" 
travels along the characteristics into (and not out of) it. When this condition is satisfied, 
a shock (discontinuity) propagates with velocity V in the Ith equation of the system. 
For I = 1, only Eq. (J2TJ is needed and A/(-) = F'(-). 

For example, the violation of the Lax condition can be used to show that tem- 
perature shock waves (under the Maxwell-Cattaneo theory) cannot propagate into a 
rigid heat conductor with temperature-dependent conductivity when the temperature 
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dependence acts to decrease the conductivity with respect to a reference value 

It remains to determine the shock velocity V. This is accomplished by returning 
to the integral form of the balance law, i.e., Eq. ([1]) and requiring it holds true for a 
shock connecting two constant states u + and u~. The result is the Rankine-Hugoniot 
"jump" conditions: 
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V(vT - u+) = F(u - ) - F(u+) (23) 
The latter set of equations can, in theory, be solved for V. We expect that there exist 
as many (not necessarily distinct) solutions for V as there are components of u. Now, 
if F(u) admits the decomposition F(u) = Au with A independent of u (as is the case 
if the nonlinear system is linearized), then Eq. f[2"3"j) reduces to 

V(u7 - u + ) = A(iT - u + ) (24) 

whence the velocity of a discontinuity in the Ith component of u is the /th eigenvalue 
of A. For weak shocks of the scalar conservation law in Eq. (fl6l) . i.e., \u~ — u + \ <C 1, 
an approximation can be made: V = \[F'(u~) + F'{u + )} + 0{\u~ — u + \ 2 ) [8, §2.6]. 
Furthermore, note that the source term R(u) does not enter the Rankine-Hugoniot 
conditions. 

For example, Christov and Jordan 14( used the Rankine-Hugoniot conditions 
to find the velocity of a heat pulse (i.e., a shock wave in the temperature) propagating 
in a rigid stationary heat conductor with weakly temperature-dependent conductivity 
under the Maxwell-Cattaneo theory. They showed that 



V/c = ^l+ e -6- (25) 

where 9 denotes a relative temperature, it has been assumed that 6 + = 0, Co is the speed 
of infinitesimal temperature waves, and e is a dimensionless parameter that determines 
the "strength" of the temperature dependence. Equation (|25i) illustrates an important 
fact about nonlinear shock wave solutions: the speed depends on the amplitude. This 
dependence is often inextricable, leading to the impossibility of deriving closed-form 
explicit expressions for both the location and magnitude of a shock wave. 
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Parabolic Equations 

Another class of PDEs relevant to problems in heat conduction and thermal stresses is 
nonlinear reaction-diffusion equations of the form 

r)9 

- = V-[D(9)Ve] + R(6) (26) 

where 9(x.,t) is the temperature change relative to some reference value, D{9) is the 
temperature-dependent diffusivity of heat, and R(9) represents either sources/sinks of 
heat in the domain or production/loss of heat due to a chemical reaction. For certain 
functional forms of D(9) and R(9), Eq. (126)) may posses solutions that behave like 
waves. 

Solution Types 

A standard technique for finding wave solutions of parabolic equations is called seeking 
a traveling wave solution (TWS). The procedure is as follows: in the one-dimensional 
context with 9 as the dependent variable, set 9(x,t) = 0(x — ct) and determine what 
and c must be. The assumption that 9 is a function of only the combined coordinate 
x — ct reduces a PDE like Eq. f l26|) to an ODE. As boundary conditions one imposes 
the growth or the value (s) of as \x\ — > oo. In multiple dimensions, the TWS ansatz 
is #(x, t) = 0(k ■ x — ct), where k is the a priori unknown direction of propagation of 
the wave. 

Though a TWS is a function of the wave coordinate x—ct, the governing equation 
lacks hyperbolicity (and, hence, finite domains of influence and dependence), so it may 
be more appropriate to call such phenomena traveling patterns rather than waves. In 
other words, parabolic PDEs may possess specific traveling patterns, however, it is 
not always possible to find solutions in the form of traveling disturbances for general 
classes of initial conditions. Yet, some (degenerate) parabolic equations admit stable 
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compact TWSs, i.e., TWSs that vanish outside of a finite spatial interval; such a TWS 



necessarily has finite domains of influence and dependence. Gilding and Kersner 
provide an advanced mathematical treatment of TWSs of nonlinear parabolic PDEs. 
A large catalog of solution types, including a variety of other solution techniques, is 
documented by |16 |. 



Temperature-dependent Conductivity 

Consider heat transfer in a rigid, stationary one-dimensional heat conductor with weak 
temperature-dependent conductivity. This leads to D(9) = 1 + e9 where |e| <C 1. 
Then, in the absence of production/loss of heat (i.e., R{9) = 0), Eq. f[2T)j) becomes 



d6_ 

dt 



d_ 

dx 



'l + et 



89 
dx 



(27) 



The TWS of the latter PDE for e < is given by Christov and Jordan |l4j : 



9{x,t) 



-l, 

w {ee^), e>e* 

where £ = x — ct, Wq(-) denotes the principal branch of the Lambert W -function 



(2? 



£* = c 1 (l + e + ln|e|), and any c > is allowed. Since the domain on which 9 is non- 
constant is semi-infinite, the TWS given by Eq. ( )28|) represents a semicompact wave 



19| with a wavefront x = £* + ct propagating to the right with phase velocity c. The 



solution given by Eq. fT28|) is illustrated in Fig. [3ta) 



Nonlinear Heat Equation 

At temperatures near absolute zero, the departure from a constant conductivity can 
be severe. It has been proposed that, as a result, the diffusivity D{9) = 9 P for met- 



als at such low temperatures 



20|. Then, for one-dimensional heat transfer in a rigid, 



stationary conductor in the absence of heat production/loss: 



d9__ d_ 
dt dx 



,09_ 

dx 
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(29) 



This PDE is sometimes referred to as the porous medium equation because it arises 
in the study of filtration of a Newtonian fluid through a porous medium; it has a 



well-developed mathematical theory [211. 
Equation $M admits the TWS 



22 



, Chap. X, §3]: 



6{x,t) 



[— pc[x — ct)] l l p , x < ct 
0, x > ct 



(30) 



for any real p, c > 0. The TWS given by Eq. Q30j) is once again a semi-compact wave 
but notice the functional form is distinctly different from that in Eq. (1281) . showing the 
case-by-case nature of wave solutions to nonlinear parabolic PDEs. If a constraint on 
the heat flux q = j_ 9(x,t) dx is imposed, then the front location is not x = ct (for 
any c > 0) but x oc 

p(p) |2^. The solution given by Eq. f)30p is illustrated in Fig. [31(b) . 
Fisher-KPP Equation 

Now, consider the case in which the diffusivity is constant, i.e., D{9) = 1, but intro- 
duce a saturating reaction term with logistic-type nonlinearity: R(9) = 9(1 —9). Equa- 
tion (126|) now becomes the so-called Fisher-Kolmogorov-Petrovskii-Piskunov (Fisher- 
KPP) equation: 



89 d 2 9 



which has the following TWS 



6(x,t) 



Jt dx 2 



23|]: 



+ 0(1- 



(31) 



[1 _|_ e (x-x -ct)/V&Y : 



5 

\/6 



(32) 



where Xq is an arbitrary constant used to set the initial front position. Unlike the 
solutions in Eqs. ( 125]) and (130|) . the TWS given by Eq. (132|) is not semicompact. Rather, 
it is a smooth wavefront, connecting two unequal constant states at x — » ±oo. Also 
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unlike the previous two examples, only one phase velocity c = 5/ a/6 is allowed under 
Eq. f )3T|) . The solution given by Eq. f )32|) is illustrated in Fig. Etc). 




Fig. 3. Traveling wave solutions 9(x,t) (superimposed at five consecutive times) of (a) the 
temperature-dependent-conductivity heat equation, (b) the nonlinear heat equation (p = 2) and (c) 
the Fisher-KPP equation. 



Further Reading 



Classical references, focusing on (amongst other things) the theory of wave solutions to 
problems in classical physics (e.g., gas dynamics, water waves, etcjjkiclude Courant 



and Friedrichs 



24( , Morse and Feshbach 



25 



Courant and Hilbert 



10] and Whitham 



26] is available on the subject, while Caviglia 



8] . An introductory textbook by Bland 

and Morro ?J cover wave propagation inhomogeneous media in detail. A review of the 



wave theory of heat conduction was written by Joseph and Preziosi 



% |5| . The topic 



has been reexamined by [27] to include recent results. Wave solutions in the study 
of thermoelasticity with finite wave speeds are discussed by Ignaczak and Ostoja- 



Starzewski 



28j. 



Modern textbooks on hyperbolic PDEs and wave solutions are necessarily fo- 
cused on the numerical as well as the theoretical aspects. LeVeque Hand Trangenstein 



ll| , who developed 



29] masterfully blend these topics. The set of lecture notes by Lax 
much of the early theory, provides a gentle introduction to the fundamental mathe- 
matical topics related to hyperbolic wave phenomena. 
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A class of wave solutions to hyperbolic equations not discussed above is those 
that are smooth in the field variables but may suffer a jump discontinuity in a higher 
derivative. Known as acceleration waves in fluid mechanics and temperature-rate waves 
in heat conduction, these "weak discontinuities" are discussed in some detail by Blanc 



26 



28 



Chap. 10] and Straughan 



27|, 



Chap. 6 §9], Ignaczak and Ostoja-Starzewski 
Chap. 4]. 

Solitons are another important category of wave solutions that were not dis- 
cussed above. They arise from a balance of nonlinearity and dispersion. In general, 
solitons are unlike the traveling wave solutions of the nonlinear parabolic PDEs dis- 
cussed above because (a) a soliton's wave speed depends on its amplitude and (b) two 
solitons can pass through each other ("interact") remaining unchanged. It has been 
suggested that soliton-like temperature waves may exist in nonlinear rigid conductors 
at low temperatures under the theory of Hetnarski and Ignaczak [30] . A detailed in- 
iroduction to the mathematical theory of solitons is presented by Ablowitz and Segur 



31], while Dauxois and Peyrard 



32] give an overview of their physics and applications. 



Cross-references 

Ordinary differential equations, Partial differential equations, Laplace and Fourier 
transforms, Special functions, Hyperbolic Heat Conduction Equation, Thermoelastic 
waves, Propagation of wavefronts in thermoelastic media, Piezothermoelasticity with 
finite wave speeds, Wave propagation in Coupled and Generalized Thermoelasticity 
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